Optical conductivity of doped Mott insulator: the interplay between correlation and 

electron-phonon interaction 
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The optical conductivity (OC) of cuprates is studied theoretically in the low density limit of 
the t-t'-J-Holstein model. By developing a limited phonon basis exact diagonalization (LPBED) 
method capable of treating the lattice of largest size 4x4 ever considered, we are able to discern fine 
features of the mid-infrared (MIR) part of the OC revealing three-peak structure. The two lowest 
peaks are observed in experiments and the highest one is tacitly resolved in moderately doped 
cuprates. Comparison of OC with the results of semianalytic approaches and detailed analysis of 
the calculated isotope effect indicate that the middle-energy MIR peak is of mostly magnetic origin 
while the lowest MIR band originates from the scattering of holes by phonons. 

PACS numbers: 71.10.Fd, 71.38.-k, 02.70.Ss, 75.50.Ee 
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The way to disclose the nature of the high temper- 
ature superconductors lies on the understanding of the 
dynamics of the holes doped into a Mott insulator [l[ . It 
is recognized that the dynamics of holes is governed by 
the interaction with magnetic subsystem as it was proved 
by the angle resolved photoemission spectroscopy of un- 
derdoped compounds [2|. There is also a growing num- 
ber of evidences that a considerable coupling to lattice 
contributes to the properties of the holes too Q . These 
major interactions are expected to leave fingerprints in 
the the OC of cuprates. However, the interpretation of 
even the basic features of the OC is controversial. 

Not to say about fine structure, there is no agreement 
on the issue of how many peaks are seen in the OC, 
both theoretically and experimentally. Initially, only the 
Drude term and MIR peak at around 0.5 eV have been 
considered as contributions coming from the dynamics 
of charged carriers 0]. Later, improved quality of the 
samples and experimental techniques gave an indication 
@, 0, B Q , and finally clearly showed [13] , that there is 
at least one more band (MIRlow band) induced by dop- 
ing in the energy range, ~ O.leV, which is just above the 
phonons energy. Moreover the analysis of the experimen- 
tal data suggests that also another contributionjpeaked 
at 1.5eV, should be also taken into account @, 0, S, 
Although there is a temptation to explain the later con- 
tribution as reminiscent of the charge-transfer peak in 
doped system, this third high energy peak (MIRhigh 
band) is observed in La2- 2; Sr 2 ;Cu04 at 1.5 eV which is 
considerably smaller than the peak energy at 2 eV, ob- 
served in undoped compound [fj. 

Even if the existence of a peak structure is recognized, 
its nature has been debated. Inability of the prototypi- 
cal t-J model, where hole moves in an antiferromagnetic 
background, to explain the experimental structure of OC 
attracted a significant interest on this problem. One pos- 



sible direction consists in considering the Hubbard model 
with moderate U. The interpretation of the mid-IR peak, 
based on a purely electronic effect (associated with the 
upper/lower-Hubbad bands and some in-gap states in- 
duced by dopi ng) . has been proposed Recently it 
was concluded [l 21 ] that the moderate U can reproduce 
the OC spectra in La2- a ;Sr :E Cu04 by U = 4eV, which 
is in sharp contrast to the t-J picture in the strong cor- 
relation limit. The other possible direction is to con- 
sider the additional electron-phonon interaction (EPI), 
i.e., t-J-Holstcin (t-J-H) model, where hole interacts also 
with dispcrsionlcss phonons. The OC of the latter model 
was calculated by several methods: exact diagonaliza- 
tion (ED) on the small \/T0 x VTO system [13j|, Self- 
Consistent Born Approximation (SCBA) with respect to 
both phonons and magnons [3| , Dynamical Mean-Field 
Theory (DMFT) for infinite dimensions [Hj], Diagram- 
matic Monte Carlo (DMC) with SCBA for magnons fiol ]. 
and ED within the Limited Functional Space (EDLFS) 
The fine structure of the OC in realistic 2D sys- 
tems can be studied only by DMC [lj| and EDLFS [jj^ 
whereas the rest of approaches encounter severe prob- 
lems [13] • In the light of this statement it is extremely 
alarming that interpretation of the low energy peak of 
OC (POClow) in and [B] is different. Magnetic ori- 
gin of POClow is concluded in [l6| while the statement 
about phononic origin in [Io| might be an error originated 
from the SCBA in magnetic channel or spin wave approx- 
imation used in [icj |. Therefore, the convincing evidence 
for the origin of the POC is an urgent and important is- 
sue towards the understanding of the basic interactions 
governing high temperature superconductors. 

In the present Letter we study theoretically the OC of 
t-t'-J-H model directly compared with the experimental 
observations supporting the vital role of EPI in cuprates. 
In addition to the inclusion of the realistic next near- 
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est neighbor hopping^' to reproduce the oberved Fermi 
surfaces in cuprates [l| , we avoid the spin wave and self- 
consistent Born approximations for the coupling to the 
spin system. By developing a LPBED method, we can 
calculate the OC of the largest ever considered 4x4 sys- 
tem. Due to the exponential growth of the basis with size 
of the system, the 4x4 lattice has considerably denser 
quantum states than the a/TO x VTO system so that it is 
possible to resolve fine structure of the OC. For the first 
time we observe, in different ranges of EPI, three peaks 
in OC, the highest one being seen because the spin- wave 
approximation is avoided. Calculating the isotope effect, 
which induces changes both in phonon frequency and ex- 
change constant J, we show that, in the weak coupling 
regime, the low energy POClow an d the middle energy 
POC are of phononic and magnetic origin, respectively. 
Furthermore, in the intermediate coupling regime the low 
energy POClow is still of purely phononic origin while 
the middle energy POC is a mixture of the lattice and 
magnetic excitations. Finally, comparison with the re- 
sults of different approximate schemes shows that the 
highest energy POChigh peak is due to incoherent tran- 
sitions into the states unaffected by lattice deformation 
associated with the hole. 

The Hamiltonian for t-t'-J-Holstein model is a sum of 
t-t'-J Hamiltonian 

H tt'j = ~t X! c i+5,c c i,° - *' X! C i+aV Ci >< 7 

i,S,cr i,8' ,a 

+ 7f^2 S i+sSi- '^^2 n i+sn t , (1) 
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EPI Hamiltonian Hh- 
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and dispersionlcss phonons Hamiltonian H j 
uiq Y^ii a l a i ( a l is the creation operator of a phonon at 
site i with frequency ojq). Here t represents the hopping 
amplitude of the site i to nearest neighbors i + 8, t' is 
the diagonal hopping amplitude to next nearest neigh- 
bors i + 6 , J is the exchange constant of the spin-spin 
interaction, a, a is the fermionic operator with excluded 
double occupancy, Si is the ^-spin operator at site i, 
and rii is the site i number operator. We introduce the 
EPI dimensionless coupling constant, A = g 2 ui /At, with 
the value A = 1 dividing the weak and strong coupling 
regimes of the Holstein model in the adiabatic limit. Be- 
low we set H = 1, t = 1, J = 0.3, if = -0.25, lu = 0.15 
and OC is in units of 2ire 2 . The one-hole ground state 
of the t-J model on 4 x 4 lattice is sixfold degenerate. 
This degeneracy between (±n/2, ±n/2), (0,7r) and (tt, 0) 
is partially removed by t' [l9| providing a four-fold degen- 
erate ground state at momentum (±7r/2, ±7r/2). Hence, 
one can naively expect that the OC should be sensitive 
to the value of t'. 

The LPBED method is based on the modified Lanc- 
zos algorithm [20j . where magnetic degrees of freedom 



are treated exactly whereas the phonon variables are ef- 
ficiently limited to a set which, as it is shown below, gives 
better results than both the Momentum Average (MA) 



approximation [2l| and SCBA. We use the translational 



symmetry associated to periodic boundary conditions, re- 
quiring that the states have a definite momentum, and 
work in the one- hole subspace with ^ S" = \ . Each 
basis vector is a linear superposition with appropriate 
phases of the 16 translational copies of a state having a 
given hole location with assigned locations of the phonon 
quanta and spin flips (hole, spin and lattice configura- 
tions are together rigidly translated). All 6435 spin con- 
figurations of 4 x 4 lattice are included. 

The real bottleneck comes from Hilbert space required 
by the phonons basis. For instance, if all the phonon 
configurations up to M=15 phonon states are included 
the size of the system is strongly limited ( \/l0 x \/10 
[HI). To circumvent this difficulty, LPBED keeps only 
two groups of phonon states. The first group are the 
lattice configurations involving only single site defor- 
mations, all the others being undeformed \phyf 1 = 

{a)) n |0) - [V^!]- 1 Ui^j |0)i- Here j = 1, ..N denotes the 
lattice sites, |0)j is the i-site phonon vacuum state, and 
all possible n = 0, 1, ...M values are limited by M = 20, 
which is shown to be enough for convergence even in the 
strong coupling regime [1 31 ] - Restriction to such basis 
makes the^method equivalent to the MA approxi- 
This is already a reasonable approxima- 
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mation 

tion, which satisfies exactly the first six spectral weight 
sum rules for the Lehmann spectral function 2l|. To 
improve the MA scheme, LPBED method includes addi- 
tional phonon states with up to three (n\ + 712 + 713 < 
3) phonons on different lattice sites Iph/^jzj'™ = 

riLX)"* [v^r 1 10),, rw, J3 io>, . such Lit is 

chosen since three-phonon basis is able to recover the 
SCBA contribution (up to three phonons) and it goes 
beyond including all other processes of the same order 
which are not present in the SCBA approach [23||. Hence, 
LPBED approach is better then MA and SCBA meth- 
ods and the only domain where it can fail is the case 
when there is a strong non-local deformation. However, 
strong deformations are realised only in the strong cou- 
pling regime where deformation in the case of local Hol- 
stein interaction is just restricted to the hole position. 
Our approach, able to improve MA, is related to that in- 
troduced in Ref . [24| , with the advantage that in LPBED 
method it is possible to calculate not only self energy of 
the quasiparticle but any correlation function. 

At T — the OC at nonzero frequency is calcu- 
lated using the Kubo expression of the OC in terms 
of the current-current correlation function a xx (uj) = 
~(Nbjy lc 3(n(uj + i?])+Il(u-iri)) with U(uj + in) = 
(tpo\jx[u + in ~H + E ]- 1 j x I Vo> • Here |Vo) is the 
ground state (k — (§,§)) with energy E , N indi- 
cates the number of lattice sites, 77 is a broadening 
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FIG. 1: (Color online) OC in (a) weak [r? = 0.025] coupling 
regime for t = —0.25 and (inset) for t = 0. OC at t' = —0.25 
(b) from intermediate to strong coupling regime for r\ — 0.1 
and (inset) 77 = 0.025. 




FIG. 2: (Color online) (a) Binding energy obtained within 
LPBED method (circles), MA method (diamonds) and SC 
scheme (squares), (b) OC in the strong coupling regime (77 = 
0.1). 

factor that shifts the poles of cr xx (u>) in the complex 
plane by replacing the 5 functions by Lorcntians, and 

jx = ietE l ^,a4+5,a^(^)^+ iet 'T, l: 8' .a C l+S' ,a C ^^'"> x - 

For A = we duly observe the well known POC in 
the t-J model at the energy 2J: the presence of t' has 
little influence on the main feature of OC of t-t'-J model 
(Fig- UK)- At weak EPI, in agreement with [n|, we detect 
POClow just above the phonon frequency. We stress 
that the existence of POClow is not related to t' hop- 
ping term, since the POClow is also observed in the t-J- 
H model (inset in Fig.[TJi). This low energy peak appears 
only at nonzero EPI and persists up to the strong cou- 
pling regime (inset in Fig.[Tb). By increasing A, the POC 
peak, which is around 2J at A = 0, shifts to higher en- 
ergies and its weight is gradually transferred to the high 
energy POChigh above u> > 2t. Hence, the EPI changes 
the spectrum of the t-J model: OC exhibits three peaks. 
The nature of these pea ks is either unknown or under 
dispute (cf. [l(| and [l6|). In the following we present 
several results unambiguously revealing genesis of these 



peaks in the weak and strong coupling regimes. 

To study the origin of the three peaks for large EPI, 
we introduce strong coupling (SC) adiabatic approach in 
which wave function is factorized into a product of nor- 
malized variational functions ^(r)) and |</>(R)) depend- 
ing on electron r and phonon R coordinates, respectively. 
The expectation value of the Hamiltonian on the state \<j>) 
provides a Hamiltonian, H e i, depending only on the elec- 
tronic degrees of freedom. It describes the t-t'-J model 
in a potential well. In Fig. [2^ we compared results of 
SC approach for ground state energies of t-t'-J model 
with data obtained by LPBED method. For A above the 
critical A c w 0.5 the results of SC, MA, and LPBED ap- 
proaches are in good agreement and, thus, the system is 
in the strong coupling regime. In the SC limit, accord- 
ing to Franck-Condon principle, the lattice is frozen in 
the ground state during the hole optical excitations and 
the OC can be calculated considering excitations of the 
hole in the static potential well formed by the lattice de- 
formation. Comparison of SC result for OC with that 
obtained by LPBED for A = 0.55 (Fig. Eb) shows that 
SC approach reproduces all three peaks. We note that in 
the SC approach both initial and final electronic states 
are calculated in the lattice potential associated with the 
ground state wave function of the hole. On the other 
hand, if one assumes that the hole in the final state re- 
leases the deformation, we obtain another curve for OC 
with only the high energy POChigh- Such approach is of- 
ten called " photoemission" (PH) process because the lat- 
tice deformation is only present in one of the counterparts 
of the states linked by current operator. Comparing re- 
sults from LPBED, SC, and PH approaches we conclude 
that POClow and POC in the strong coupling regime 
represents hole transitions between states within the self- 
consistent potential well generated by the phonons. To 
the contrary, POChigh is associated with transitions into 
states which are intact by the EPI driven lattice deforma- 
tion. Here we note that the association of the theoretical 
POChigh with experimental MIRhigh band should be 
done with care, since both structures are located at en- 
ergies which are near the limit of applicability of the t-J 
model, where one is required to use methods which can 
handle the initial unreduced three-band Hubbard model 
[25| . We stress that the above discussed scenario is valid 
also in the more simple t-J-Holstein model. 

Another way to establish the nature of the peaks of 
OC is to study changes of OC with small variations of the 
phonon frequency uiq and/or exchange constant J. Such 
study does not serve solely as Gedanken experiment but 
establishes how the spectra will be altered by the isotope 
substitution (IS). First of all, changes of the oxygen from 
16 to 18 induces modifications in the values of coupling 
constant g and phonon frequency cJo- Phonon frequency 
ljq = ^/k/M is expressed in terms of k, which is the 
restoring force for length unit of the local oscillators, and 
M, which is the mass of oxygen atoms surrounding the 
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FIG. 3: (Color online) Effect of the O IS on OC in the weak 
[(a), (bi), and (b2)], and the intermediate coupling regime 
[(c), (di), and (d 2 )]. 



creasing ui [28|. Thus, the decrease of lu increases the 
energy of POC. To the contrary, decrease of J tends to 
soften the POC. Hence, joint influence of both effects 
leaves the peak at the same position. To confirm this 
conclusion we repeated the calculations setting to zero 
the change of J and found (not shown) the increase of 
POC energy. Hence, we conclude that in the intermedi- 
ate coupling regime the middle POC is of mixed origin. 

In conclusion, we developed a novel method capable 
of studying the fine structure of the OC of the t-t'-J- 
Holstcin model and found that the influence of the diag- 
onal hopping t' on OC is surprisingly little. For the first 
time we revealed a 3-peak structure and established the 
origin of these peaks in the whole range of hole-phonon 
couplings. We also predicted the influence of the IS on 
the OC which can be compared with experiment after the 
problem of residual 16 O in the matrix of 18 O, strongly in- 
fluencing low energy part of OC will be solved. 

A.S.M. was supported by RFBR 07-02-00067a, N.N. 
by Grant-in-Aids No. 15104006, No. 16076205, No. 
17105002, No. 19048015, and NAREGI Japan; G.D.F., 
V.C. and C.A.P. received financial support from Research 
Program MIUR-PRIN 2007. 



Cu ion in the Cu02 plain. IS changes the values of wo and 
jtow * = ljot/M/M* and g* = g(M*/M) 1 / 4 with value 
of A independent on isotope. In particular, the relative 
shift of ojq is about 6%. The second possible effect of IS is 
the decrease of the antifcrromagnctic exchan ge c onstant 
J in compounds with the apical oxygen 26|, 27 1, driven 
by its vibrations out of plane. With the IS the value of J 
is reduced by about 1% and in the following we assume 
AJ/J = -0.01. 

In Fig. 3 we present the changes of OC induced both 
by change of the in-plane oxygen mass and exchange 
constant J. The contribution of these changes are well 
distinct with only one exception. In the weak coupling 
regime (Fig. [3^ and Fig. [Hp! ) POC at frequencies around 
2 J is shifted down by about 1% indicating the clear mag- 
netic origin of this peak in the weak coupling regime. 
To the contrary, POClow, with frequency around Wo, is 
shifted down about 6% (Fig. \3lp2), that is just softening 
of loq induced by IS. Thus, we get one more confirma- 
tion of the phononic origin of theoretical POClow and 
experimental MIRlow band, in agreement with Ref. [loj ] . 

One gets the same conclusion about POClow from 
the IS in the intermediate coupling regime (Fig. and 
Fig. Hi). Indeed the POClow is shifted down by 6% 
again. On the contrary, behavior of the middle energy 
POC differs from that in the weak coupling regime since 
it reduces its intensity and almost does not move with 
the IS. It is known that the middle POC energy increases 
by increasing A/A c [lj| • Moreover, A c decreases by dc- 
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